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ABSTRACT 

New spectroscopic observations of 36 H II regions in NGC 4258 obtained with the Gemini telescope are com- 
bined with existing data from the literature to measure the radial oxygen abundance gradient in this galaxy. The 
[O III] A4363 auroral line was detected in four of the outermost targets (17 to 22 kpc from the galaxy center), 
allowing a determination of the electron temperature T e of the ionized gas. From the use of different calibra- 
tions of the R23 abundance indicator an oxygen abundance gradient of approximately -0.012±0.002 dexkpc" 1 
is derived. Such a shallow gradient, combined with the difference in the distance moduli measured from the 
Cepheid Period-Luminosity relation by Macri et al. between two distinct fields in NGC 4258, would yield 
an unrealistically strong effect of metallicity on the Cepheid distances. This strengthens the suggestion that 
systematic biases might affect the Cepheid distance of the outer field. Evidence for a similar effect in the 
differential study of M33 by Scowcroft et al. is presented. A revision of the transformation between strong- 
line and T^-based abundances in Cepheid-host galaxies is discussed. In the T e abundance scale, the oxygen 
abundance of the inner field of NGC 4258 is found to be comparable with the LMC value. 

Subject headings: galaxies: abundances — galaxies: ISM — galaxies: individual (NGC 4258) 



1. INTRODUCTION 

The geo metric distance to the spiral galaxy NGC 4258 ob- 
tained by |Herrnstein et al.| ( [T999] 7.2 ± 0.5 Mpc) from the ke- 
plerian motion of water masers in its nucleus represents a cor- 
nerstone of the extragalactic distance ladder. The prospect 
of reducing the uncertain ty in the distance from 7% to 3% 
(Humphre ys et al.| |2008) opens the possibility of adopting 
NGC 4258 a s a new anchor of the extragalactic distance scale 
( |Riess et"al . 2009a), bypassing the use of the LMC, whose 
distance remains uncertain at the 5-10% level. Moreover, the 
availability of a geometric distance for a galaxy located well 
beyond the Local Group offers the opportunity to cross-check 
the results obtained from primary and secondary distance in- 
dicators at distances that are comparable to those of the cal- 
ibrators of far-reaching standard candles, such as the type la 
supernova peak brightness. Distances to NGC 4258 that are 
consistent with the geometrically determined one have been 
obtained from the Ce pheid Period-Luminosity (P-L) relation 
( Newman et al.||2001[) and the Tip of the Red Giant Branch 
(TRGB, |Mager et al.|2008[ ). From the Cepheids observed with 
the Hubble Space Telescope (HST) by Macri et aL] ( |2006[ l var- 
ious authors, using different treatments of the extinction and 
metallicity corrections, together with the adoption of differ- 
ent P-L relations, have obtained distance moduli that differ at 
the 0. 1 mag level, but that are still compatibl e, within the re- 
ported uncertainties, with the m aser dist ance (|An et al.|2007[ 
|van Leeuwen et al.|2007l|Benedict et al.|2007l|Turner|2 010). 

The effect of metallicity on the Cepheid P-L relation 
is a crucial but still debated issue. A persistent discrep- 
ancy between theoretical predictions from pulsation models 

1 Based on observations obtained at the Gemini Observatory, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under a cooperative agreement with the NSF on behalf of the Gem- 
ini partnership: the National Science Foundation (United States), the Sci- 
ence and Technology Facilities Council (United Kingdom), the National 
Research Council (Canada), CONICYT (Chile), the Australian Research 
Council (Australia), Ministerio da Ciencia e Tecnologia (Brazil) and Min- 
isterio de Ciencia, Tecnologia e Innovation Productiva (Argentina). 



( Fiorentino et al.p007||Bono et al.|2008jl and em pirical deter- 
minations ( |Macri et al.| |2006; Scowcroft et al. 2009) lingers 
on. In the HST Key Project ( [Freedman et al.||2001| ) a metal- 
licity correction factor 7 = (5/io/<5 logZ = -0.20 mag dex" 1 was 
applied to the galaxy distance moduli, in the sense of making 
them larger for metallicities exceeding the LMC value. The 
estimated contribution of the metallicity to the pres ent 5% un- 
certainty on the H ubble constant Hq is about 2.5% (Freedman 
& Madore 2010). Current efforts striving for an accuracy bet- 
ter than 5% on Hq attempt to circumve nt this issue by sele ct- 
ing galaxies with similar metal content ( |Riess et al.|2009b I. 

Given the importance of NGC 4258 in current investiga- 
tions of the extragalactic distance scale, in this paper the oxy- 
gen abundance of its interstellar medium (a proxy for the 
Cepheid metallicity) is revisited. One of the initial goals of 
this new study was to determine oxygen abundances from 
direct electron temperature (T e ) determinations in Hll re- 
gions. This requires deep integrations with large aperture 
telescopes, in order to detect the weak [O III] A4363 auroral 
line, which provides, in combination with the much stronger 
[O III] AA4959, 5007 nebular lines, the crucial temperature di- 
agnostic for ionized nebulae. Despite the fact that NGC 4258 
has been the subject of previous abundance studies based on 
its H II region population, the A4363 line has remained un- 
observed until the present time. Secondly, it was felt that an 
investigation based on new spectroscopic data was warranted 
in order to verify the recent determina tion of the metallicity 
effect on the P-L relation obtained by Macri et al. (2006J from 
a differential study of Cepheids located in two fields at differ- 
ent distances from the galaxy center (and therefore different 
metal content). This study adopted metallicities deter mined 
from the calibration of the R23 strong-line indicator by |Zarit-| 



I sky et aL] ( fl994[ ), applied to H II region emission-line dat a 
obtained by the same authors and Oey & Kennicutt (1993). 
However, a recent reanalysis of these data by |Bono et al.| 



( 2008 ) suggested that the abundance gradient might need to 
be revised. Lastly, during the course of this work some incon- 
sistencies were found in the r e -based abundance scale adopted 
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in several published studies of extragalactic Cepheids (Sakai 
|etal.|2 004, Sah a et al.|2006| l, and these will also be addressed 
for the benefit of future studies. 

It is important to keep in mind that the determination of ex- 
tragalactic chemical abundances from nebular spectroscopy is 
still afflicted by large systematic uncertainties, despite the ma- 
turity of the field (see Bresolin|2008 for a recent review). A 
well-established discrepancy exists between abundances ob- 
tained from the classical, direct r e -based method ( |Menzel| 
|et al.|[T94T| and the theoretical p redictions of photoioniza- 
tion model grids (McGaugh 1991). In addition, the study of 
metal recombination lines provides oxygen abundances that 
lie between the direct (from collisionally excite d lines) and 
the theoretical values (Garcfa-Rojas et al. 2007). Generally, 



Table 1 

NGC4258: adopted parameters 



Parameter 



Value 



it is believed that differential analysis methods are unaffected 
by these systematic differences. However, in the case of the T e 
vs. strong-line methods in the high metallicity regime (around 
the solar ox ygen abundance 12 + log(O/H) = 8.69, Asplund 
|et al.|20 09) this conclusion is still based on a limited number 
of experimental data points. 

2. OBSERVATIONS AND DATA REDUCTION 

New spectra of H II regions in NGC 4258 were ob tained 
with the G emini Multi -Object Spectrograph (GMOS, |Hook 
et al.|2004 1 at the Gemini North facility on Mauna Kea. The 
targets were selected from Ha narrow-band images of two 
5(5x5^5 GMOS fields, centered approximately 2(3 NW and 
6(5 SE of the center of the galaxy, obtained between January 
and March 2010. The spectroscopic data were acquired in 
queue mode between March and June 2010, using two multi- 
object masks, one per field, with l('5-wide slitlets. The see- 
ing conditions were around 0"8 during the observations of 
the NW field, and around 0"5 for the SE field. In order to 
minimize the effects of the atmospheric differential refraction 
the data were acquired at airmasses smaller than 1.20. Three 
3000 s exposures were secured for each of the two fields us- 
ing the B600 grating, which provided spectra covering at least 
the 3500-5100 A wavelength range at a spectral resolution of 
^5.5 A. For some of the targets, depending on their spatial 
distribution, the spectral coverage extended up to ^5900 A. 

iRAtj^] routines contained in the gemini/gmos package 
were used for cosmic ray rejection, electronic bias subtrac- 
tion, flat field correction and wavelength calibration of the raw 
data frames. Observations of the spectrophotometric standard 
Feige 34 yielded the flux calibration. Thanks to the similarity 
in airmass and atmospheric conditions between the three dif- 
ferent frames obtained for each individual target, the emission 
line intensities of the single extracted spectra are very consis- 
tent with each other. The working version of the spectra was 
then obtained by averaging the three spectra corresponding to 
each individual slit. 

The final H II region sample comprises 36 objects, equally 
divided between the NW and the SE fields. Their celestial co- 
ordinates are summarized in Table [2] (where objects are listed 
in order of decreasing declination), together with their galac- 
tocentric distances in kpc, deprojected using the geometric 
parameters in Table [T] Fig. [T] shows the location of the spec- 
troscopic targets on an Ha image of NGC 4258 (circles num- 
bered 1 to 36). It can be seen that all but three of the objects 

2 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 



R.A. (J2000.0) 12 18 57.50" 

Decl. (J2000.0) 47 18 14.3" 

Morphological type SAB(s)bc b 

Distance 7.2 Mpc c 

#25 9f31 b (19.50 kpc) 

Inclination 72° d 

Position angle of major axis 150° d 

I!'; 8.53 b 

M° -20.76 

References. — (a) Herrnstein et al. (2005); (b) de Vaucouleurs 
|eTaT]pWl) ; (c) |Herrn'tein el :al.Hl^ ' T fifijn AlbadaHl^Oj 

lie within the R25 radius (i.e. within the ellipse in the figure). 
The H II regions observed in the NW field sample a portion 
of the northern inner spiral arm of the galaxy, together with 
a few objects close to the galaxy center and in the northern 
outer arm. The H II regions in the SE field are mostly located 
in the southern outer arm (a description of the complex spiral 
morphology of NGC 4258 is given by |Courtes et al.|1993| . 

2.1. Line fluxes and comparison with previous work 

The emission line intensities were measured with the 
splot program in IRAF by integrating the fluxes un- 
der the line profiles. The metal lines [O II] A3727, 
[O in] AA4959, 5007 and the H/3, H7 Balmer lines were 
present in the spectra of all the 36 targets in Table [2] The 
weaker [Ne III] A3868 and H<5 lines were measured inl9 and 
33 H II regions, respectively. The line intensities were cor- 
rected for interstellar reddening by assuming a case B in trinsic 
H7/H/3 ratio of 0.47 at T e = 10 4 K and the |Seaton| < fl979| l red- 
dening law. The resulting reddening-corrected emission lines, 
normalized to Hf3 = 100, are presented in Table|3| The line er- 
rors reported in the table account for uncertainties in the flat 
fielding, the flux calibration, the positioning of the continuum 
level around each line and the extinction coefficient c(H/3). 

A few of the brightest objects in the sample have already 
been studied spectroscopically by other authors, in particu - 
lar six targets are in c omm on withfOey & Kennicutt ( |1993[ ), 
|Zaritsky et~aL] ( |1994| i and |Bresolin e t al. l (|1 9 99|l,~and these 
are identified in column 5 of Table 121 Fig. [2] compares the 
[0 11] A3727 (circles) and [O III] AA4959, 50071 (triangles) line 
fluxes measured in the present study with those measured by 
these other authors. The agreement is found to be excellent: 
the mean ratio between the two sets is 1 .01 ± 0.04. 



3. OXYGEN ABUNDANCES 

3.1. Electron temperatures and direct abundances 

One of the goals of this new spectroscopic investigation of 
NGC 4258 was to measure, for as many H II regions as pos- 
sible, the electron temperature of the nebular gas via the stan- 
dard method that involves the detection of the [O III] A4363 
auroral line. No previous detection of this line has been re- 
ported in the literature for NGC 4258. In this new study, the 
A4363 line was measured for four H II regions (numbers 24, 
30, 34 and 36 in Table|2]i, all of them located at large galacto- 
centric distances and along the southern outer spiral arm. The 
nebular package in IRAF was used to derive electron tem- 
peratures (adopting an ele ctron density of 10 2 c m~ 3 and the 
atomic parameters used by Bresolin et al.|2009" 1, and the + , 
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Figure 1. Location of the H II regions studied in this work on a narrow-band Ha image of NGC 4258 (courtesy L. van Zee). The H II regions are numbered in 
order of decreasing declination. Letters identify objects drawn from the literature. The ellipse represents the location of the projected 7?25 radius (R25 = 9'3 1 ). 
The inner and outer parts of the galaxy are shown with different stretchings, in order to facilitate the identification of the targets. 



O" 1 " 1 " ionic abundances. The resulting temperatures and total 
O/H abundance ratios are presented in Table |4] The elec- 
tron temperatures and the 12 + log(0/H) abundances are re- 
markably similar among the four targets, with mean values of 
1 1400 K and 8.20, respectively. 

3.2. Strong-line abundances 

For the majority of the H II regions in our sample we need to 
resort to the use of strong-line abundance indicators, since the 
auroral lines remain undetected. The difficulties and poten- 
tial pitfalls of the se methods have been described at length in 
the literature (see Bresoli n et al.|2009| for a recent overview) 
and will not be repeated here. In brief, with these methods 
the value of the oxygen abundance is related to the strength 
of metal lines in the optical nebular spectra either via an em- 
pirical calibration obtained from H II regions with auroral line 
detections ('empirical' method) or via the results of theoret- 
ical photoionization models. The abundances derived from 
the different procedures and calibrations vary in a systematic 
way by large factors, with empirical abundances typically oc- 
cupying the low end of the range. While this situation is cur- 



rently unexplained, relative abundances are generally consid- 
ered to be rather robust. Investigations of additional present- 
day abundance indica tors in galaxies (e.g. blue supergiants, 
Kudritzki et al. 2008) are presently being carried out, with 
the goal of establishing the reliability of the absolute nebu- 
lar abundances obtained by the different calibrations. Recent 
work in the nearby galaxies NGC 300 and M33 has shown 
that a good agreement is obtained between the supergiant and 
H II region abundances when the latter are based on the direct, 
[O III] A4363-based method ( jBresolin et al.|2009||20T0"l l. For 
this reason, it is important to be able to anchor the strong-line 
abundances obtained in more distant galaxies to the values 
provided by the detection of the [O III] A4363 auroral line, as 
done in the present investigation of NGC 4258. While rela- 
tive abundance determinations seem to be, to first order, unaf- 
fected by the choice of tec hniq ue used to measure them (but 
see last paragraph of Sec. [67Tb, it is also true that it is im- 
portant to establish the absolute nebular abundances in order 
to constrain the chemical evolutionary status of star-forming 
galaxies, a fact that is often overlooked. 
For this study of NGC 4258 it was decided to investi- 
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Table 2 

Observed H II region sample 



ID 


R.A. 


Decl. 


R 


other ID 




(J2000.0) 


(J2000.0) 


(kpc) 




(1) 


(2) 


(3) 


(4) 


(5) 



1 . . 

2.. 
3.. 
4.. 
5.. 
6. . 
7.. 
8.. 
9.. 
10. 
1 1 . 
12. 
13. 
14. 
15. 
16. 
17. 
18. 
19. 
20. 
21 . 
22. 
23. 
24. 
25. 
26. 
27. 
28. 
29. 
30. 
31. 
32. 
33. 
34. 
35. 
36. 



12 18 53.31 
12 18 55.26 
12 18 55.69 
12 18 48.52 
12 18 49.39 
12 18 51.05 
12 18 54.13 
12 18 55.20 
12 18 55.94 
12 18 56.55 
12 18 57.37 
12 18 58.93 
12 18 49.51 
12 18 58.52 
12 18 55.87 
12 18 56.19 
12 18 56.28 
12 18 43.86 
12 19 03.56 
12 19 05.18 
12 18 59.29 
12 18 56.96 
12 19 03.01 
12 18 59.09 
12 19 23.88 
12 19 24.87 
12 19 07.63 
12 19 09.78 
12 19 13.34 
12 19 02.38 
12 19 12.68 
12 19 17.94 
12 19 21.18 
12 19 19.65 
12 19 11.33 
12 19 19.49 



47 22 35.5 
47 22 09.5 
47 21 54.3 
47 21 24.8 
47 21 08.2 
47 20 53.3 
47 20 50.3 
47 20 35.5 
47 20 25.5 
47 20 15.2 
47 20 04.3 
47 19 30.6 
47 19 20.4 
47 19 05.0 
47 18 18.7 
47 18 11.2 
47 17 46.0 
47 17 25.7 
47 15 10.6 
47 15 00.8 
47 14 14.8 
47 13 43.3 
47 13 37.6 
47 13 33.8 
47 13 03.2 
47 12 59.6 
47 12 40.2 
47 12 16.0 
47 12 01.6 
47 1 1 40.7 
47 11 39.6 
47 1 1 22.0 
47 11 01.0 
47 10 50.6 
47 10 46.4 
47 09 56.0 



13.66 
13.36 
12.70 
7.58 
6.92 
6.48 
7.62 
7.36 
7.23 
7.01 



34C, (-032+141), B3 
39C, (-025 + 130), B4 
44C, (-019+121), B5 



6.94 


52C, (-010+110) 


6.09 


61C 


5.43 


18C 


4.10 


56C 


1.43 


38C 


1.48 


42C 


2.89 


43C 


16.35 


3C 


7.93 


82Ca, (+039-185) 


7.93 


83Ca 


13.97 


7Sa 


17.79 


ISa 


13.79 




16.74 




16.55 


29S 


17.21 


32S 


14.79 


2S 


15.28 


4S 


15.04 




21.56 


ZKH 8 


16.30 




16.36 




17.33 




17.61 


16S 


19.72 




19.97 





Note. — Units of right ascension are hours, minutes and seconds, and units of 
declination are degrees, arcminutes and arcseconds. Col. (1): Hll region identi- 
fication. Col. (2): Right ascension. Col. (3): Declination. Col. (4): Deprojected 
galactocentric distance in kpc. Col. (5): Main identification from Courtes et a k|(1993), 
with additional objects from Zaritsky et al. 1 1994 =ZKH), Oey & Kennicutt 1 1993 
coordinates in brackets) and |Bresolin et al.»'lW|=B). 



gate the abundances obtained from the strong-li ne indicator 
R23 = ([Oh]A3727 + [Offl]AA4959,5007)/H/3 ( |Pagel et al. 
1979|l. Two dif ferent calibrations are adopted: t he one "5y 
McGaugh] (j 1 99 1 1 in the analytical form given by |Kuzio de 
Naray et al. 2004 ), which is based on a grid of photoioniza- 
tion models, and the one by [Pilyugin & Thuan ( 2005 1. The 
latter was selected as an alternative because it is essentially 
calibrated from auroral Unerased abundances, and thus it 
should provide results that are in agreement with those ob- 
tained from the [O III] A4363 detections. Both calibrations 
account for the excitation of the ionized gas, via the pa- 
rameters y = log(03/02) and P = 03/(03 + 02), respectively, 
where 03 = [O in] AA4959, 5007 and 02 = [O n] A 3727. 

The double-valued nature of R23 represents a major source 
of uncertainty when it is difficult to ascertain whether the 
targets fall onto the upper (high O/H) or lower (low O/H) 
branch of the indicator. This issue is commonly resolved 
by considering emission line ratios that are mono tonic with 
O/H, such as [Nll]/Ha and [Nll]/[Oll]. Our spectra do not 
cover the region around Ha, and therefore we can neither 
use indicators based on the strength of the [N II] A6583 line 
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Figure 2. Comparison between the [O II] A3727 and [O III] AA4959, 5007 
reddening-corrected fluxes measured in this work and in previous publica- 
tions. The dashed line represents the one-to-one line. 



(e.g. N2 = log([Nll]/Ha, Pettini & Pagel 2004), nor use this 
information to select the appropriate R23 branch. Neverthe- 
less, it is possible to select the upper branch calibration for 
the full sample with high confidence, based on the following 
considerations: 

(a) Fig. [3] shows the relation between R23 and O/H for Hll 
regions in a number of spiral and irregular galaxies, for 
which the oxygen abundance was derived from high-quality 
measurements of the [O III] A4363 line and extracted from 
the following references: 



NGC 300: 
M101: 
NGC 2403: 

M33: 

IC 1613: 
H II galaxies: 



Kennicutt et al. K003 i 

Garnett et al.|ifT997ll 
. 1 11 1^1 



bsteban et ; 



Bresolin et al. (2009) 



If l |Z009j) 



Bresolin et al. tZUT t 
Fsteban et al 2TJ09| 



Bresolin et 



izotov et al. ( 1994 



Guseva et al 



al ( 



2006) 



2009 



Fsteban et al |pDQg] > 



M31,NGC 2363, NGC 1741 
NGC 4395, NGC 4861 

In selecting the spiral galaxy sample preference was given 
to investigations carried out by the author and his collabora- 
tors, in order to ensure consistency in the reduction and data 
analysis procedures. The data plotted in Fig. [3] clearly de- 
lineate the two R23 branches, with the turnover point located 
around 12 + log(O/H) = 8.0. Most of the Hll galaxies shown 
here fall below the turnover point and lie on the lower branch. 
Virtually all of the H II regions in spiral galaxies are located 
on the upper branch (the notable exception is represented by 
region SDH 323 in M101, which lies at a very large distance 
from the galaxy center). The lines drawn in Fig. [5] are poly- 
nomial fits to the data, for 12 + log(0/H) above anabelow the 
turnover point, and are shown to guide the eye in identifying 
the two R23 branches. The four squares in Fig. [3] represent 
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Table 3 

Reddening-corrected line fluxes and strong-line oxygen abundances 



ID 


[On] 


[Ne III] 


Hi 


[Om] 


[Om] 


log(R 23 ) 


F(H/3) c(H/3) 


12 + log(0/H) 




3727 


3868 


4101 


4959 


5007 




(erg s~' cnT 2 ) (mag) 


M91 P05 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) (9) 


(10) (11) 



10. 

11 . 

12. 
13. 
14. 
15. 
16. 
17. 
18. 
19. 
20. 
21 . 
22. 
23. 
24. 
25. 
26. 
27. 
28. 
29. 
30. 
31. 
32. 
33. 
34. 
35. 
36. 



273 
263 
290 
248 
249 
247 
219 
256 
182 
205 
216 
189 
206 
234 
199 
243 
222 
246 
270 
329 
194 
274 
422 
213 
113 
374 
355 
286 
337 
149 
314 
287 
241 
314 
326 
193 



± 28 
± 27 
±29 
±25 
±25 
±25 
± 22 
±26 
± 18 
±21 
± 22 
± 19 
±21 
± 24 
±20 
± 24 
± 22 
±25 
± 27 
± 33 
±20 
± 28 
±44 
± 22 
± 13 
±39 
±36 
±31 
±35 
± 15 
±32 
±30 
± 27 
±32 
±33 
±20 



13.9 ± 1.8 
12.5 ± 1.7 
10.5 ± 1.6 



2.9 ± 0.4 

3.6 ±0.4 
2.8 ±0.5 

3.7 ±0.9 
6.7 ± 1.7 
6.1 ± 0.9 



7.3 ± 1.0 
6.5 ± 0.8 
4.9 ±0.8 
8.7 ± 1.1 
23.6 ± 2.3 
13.9 ±2.5 



47.4 ± 4.5 



; ± 1.9 



23 ± 2 
26 ±2 
26 ±2 
26 ±2 

25 ± 2 
20 ±2 

26 ±2 

29 ±2 

24 ±2 

24 ±3 

25 ±2 
24 ±2 

27 ± 2 

24 ±2 

25 ± 2 

26 ±2 
26 ±2 

28 ±2 
26 ±2 

30 ±2 
26 ±3 
26 ±2 

24 ±4 

25 ±2 
23 ±3 



27 ± 3 

25 ± 4 
30 ± 4 

47.2 ± 4.5 29 ± 2 

26 ± 3 



26 ±2 
23 ± 3 

32.9 ± 4.3 25 ± 2 



63 ±4 


190 ± 11 


0.72 ± 0.04 


3.3 


Xl0~ 


5 


0.23 


8.67 ± 0.04 


8.35 ± 0.05 


60 ±4 


178 ± 11 


0.70 ± 0.04 


3.2 


xl0~ 


5 


0.18 


8.69 ± 0.04 


8.37 ± 0.05 


58 ± 3 


166 ± 10 


0.71 ±0.04 


3.3 


xl0~ 


15 


0.12 


8.67 ± 0.04 


8.33 ± 0.05 


44 ± 3 


130 ± 8 


0.62 ± 0.04 


1.7 


xl0~ 


15 


0.28 


8.75 ± 0.03 


8.40 ± 0.05 


28 ± 2 


80 ± 5 


0.55 ± 0.04 


4.3 


xl0~ 


15 


0.42 


8.81 ±0.03 


8.36 ± 0.05 


21 ± 1 


65 ±4 


0.52 ± 0.04 


2.6 


xl0~ 


15 


0.17 


8.82 ± 0.03 


8.34 ± 0.05 


14 ± 1 


42 ± 3 


0.44 ± 0.04 


2.1 


xl0~ 


15 


0.33 


8.88 ± 0.03 


8.34 ± 0.04 


22 ± 1 


67 ±4 


0.54 ± 0.04 


2.6 


xlO- 


4 


0.41 


8.81 ±0.03 


8.33 ± 0.05 


38 ± 2 


115 ± 7 


0.53 ± 0.04 


5.5 


xl0~ 


4 


0.27 


8.84 ± 0.03 


8.53 ± 0.04 


30 ±2 


91 ± 5 


0.51 ±0.04 


1.7 


xl0~ 


4 


0.50 


8.84 ± 0.03 


8.46 ± 0.04 


23 ± 1 


68 ±4 


0.49 ± 0.04 


1.3 


xl0~ 


14 


0.20 


8.86 ± 0.03 


8.41 ± 0.04 


69 ±4 


205 ± 12 


0.67 ± 0.03 


4.7 


xl0~ 


15 


0.44 


8.74 ± 0.03 


8.49 ± 0.04 


40 ±2 


119 ± 7 


0.56 ± 0.04 


7.8 


xl0~ 


15 


0.05 


8.81 ± 0.03 


8.48 ± 0.04 


29 ± 2 


87 ± 5 


0.54 ± 0.04 


6.1 


xl0~ 


15 


0.35 


8.81 ±0.03 


8.40 ± 0.04 


24 ± 1 


72 ±4 


0.47 ± 0.04 


5.1 


xl0~ 


5 


0.00 


8.87 ± 0.03 


8.46 ± 0.04 


50 ±3 


148 ± 9 


0.64 ± 0.04 


7.9 


xl0~ 


15 


0.37 


8.74 ± 0.03 


8.41 ± 0.05 


44 ±3 


128 ± 8 


0.59 ± 0.04 


6.5 


xlO- 


15 


0.16 


8.78 ± 0.03 


8.45 ± 0.04 


38 ±2 


112 ± 7 


0.60 ± 0.04 


4.2 


xl0~ 


15 


0.00 


8.78 ± 0.03 


8.40 ± 0.05 


43 ± 2 


129 ± 8 


0.64 ± 0.04 


5.5 


xl0~ 


15 


0.39 


8.73 ± 0.04 


8.36 ± 0.05 


74 ±4 


223 ± 13 


0.80 ± 0.04 


9.0 


xl0~ 


15 


0.90 


8.58 ± 0.05 


8.25 ± 0.06 


56 ±3 


176 ± 10 


0.63 ± 0.03 


1.9 


xlO~ 


5 


0.66 


8.76 ± 0.03 


8.50 ± 0.04 


26 ±2 


74 ± 4 


0.57 ± 0.04 


2.1 


xlO~ 


15 


0.20 


8.78 ± 0.04 


8.30 ± 0.05 


48 ±3 


142 ± 8 


0.79 ± 0.04 


1.0 


xlO~ 


15 


1.14 


8.56 ± 0.05 


8.09 ± 0.07 


157 ± 9 


465 ± 28 


0.92 ± 0.03 


1.6 


xlO~ 


15 


0.26 


8.50 ± 0.04 


8.29 ± 0.05 


104 ±6 


303 ± 18 


0.71 ±0.03 


6.0 


xlO~ 


16 


0.35 


8.73 ± 0.03 


8.54 ± 0.03 


29 ±2 


91 ± 5 


0.70 ± 0.04 


8.2 


xlO- 


16 


0.35 


8.65 ± 0.05 


8.13 ±0.06 


19 ± 1 


58 ± 3 


0.64 ± 0.04 


1.5 


xlO~ 


15 


0.43 


8.71 ±0.05 


8.11 ± 0.06 


57 ±4 


171 ± 10 


0.71 ±0.04 


3.9 


xlO~ 


16 


0.12 


8.67 ± 0.04 


8.33 ± 0.06 


56 ±3 


167 ± 10 


0.75 ± 0.04 


6.2 


xlO~ 


<i 


0.63 


8.62 ± 0.05 


8.24 ± 0.06 


169 ± 10 


501 ± 30 


0.91 ± 0.03 


1.3 


xlO~ 


4 


0.28 


8.54 ± 0.04 


8.34 ± 0.04 


76 ±4 


231 ± 14 


0.79 ± 0.04 


7.8 


xlO~ 


16 


0.08 


8.59 ± 0.04 


8.28 ± 0.06 


14 ± 1 


42 ± 3 


0.54 ± 0.04 


8.5 


xlO~ 


16 


0.07 


8.80 ± 0.04 


8.20 ± 0.05 


118 ± 7 


353 ± 21 


0.85 ± 0.03 


3.2 


xlO- 


16 


0.08 


8.56 ± 0.04 


8.33 ± 0.05 


80 ±5 


237 ± 14 


0.80 ± 0.03 


5.5 


xlO~ 


15 


0.20 


8.58 ± 0.04 


8.28 ± 0.06 


52 ±3 


157 ± 9 


0.73 ± 0.04 


8.0 


xlO" 


16 


0.00 


8.64 ± 0.04 


8.26 ± 0.06 


169 ± 10 


492 ± 29 


0.93 ± 0.03 


6.0 


xlO" 


16 


0.00 


8.50 ± 0.04 


8.29 ± 0.05 



Note. — The line fluxes are in units of H/3 = 100. F(H/3) in column (8) is the measured H/3 flux, uncollected for extinction. The oxygen abundances in columns (10) and 
(11) are obtained with the R23 calibrations by McGaugh 1 1991] =M91) andfPilyugin & Thuan ( 2005 =P05), respectively. 



Table 4 

Direct method abundances 

ID [OlII]A4363 T e (K) 12 + log(0/H) 
(1) (2) (3) (4) 

24 4.5 ± 0.7 11350 ± 850 8.23 ± 0.10 

30 5.0 ±0.4 11490 ± 530 8.19 ±0.06 

34 1.9 ± 0.2 10640 ± 510 8.23 ± 0.07 

36 5.6 ± 1.0 12030 ± 1050 8.16 ±0.11 

Note. — The reddening-corrected intensity of [O III] A4363 in column 2 is 
in units of Hfi = 100. 



the H II regions in NGC 4258 for which a direct abundance 
was obtained (see Table [4]). It can be seen that these objects 
are placed on the upper branch. Considering the fact that 
these H II regions are among the most distant from the cen- 
ter of NGC 4258, and that galactic abundance gradients are 
always measured to be monotonic, with abundances increas- 
ing toward the central regions, this result strongly suggests 
that the remaining Hll regions in this galaxy (which have 



smaller R23 values, as shown in column 7 of Table [3} have 
larger abundances, and should therefore also be placed on the 
upper branch. 

(b) all of the Hll regions without [Olll]A4363 detec- 
tions have log(R23) < 0.80, i.e. significantly lower than the 
turnaround point where the two branches intersect. From 
Fig.[3]we can infer that misplacing a lower branch object onto 
the upper branch would increase the estimated O/H ratio by 
approximately 0.5 dex or more. This would easily show up 
in a plot of the galactocentric abundance gradient, under the 
assumption that the gradient is monotonic and without abrupt 
discontinuities. If we assume upper branch objects only for 
our sample, the gradient in NGC 4258 behaves 'normally'. If 
instead we assume that some of the objects at high R23 values 
belong to the lower branch, we obtain very large discontinu- 
ities in the abundance gradient, which can be reasonably ex- 
cluded for H II regions located within the optical boundaries 
of normal galaxies. 

(c) For four of the targets (4, 5, 13 and 34 in Table [3]) the 
auroral [N II] A5755 line was within the covered wavelength 
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Table 5 

Additional H II regions from the literature 
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Figure 3. Log(R23) vs. 12 + log(0/H) for a sample of extragalactic H II re- 
gions drawn from the literature. The oxygen abundances were all derived 
using the direct method, i.e. using the [O III] A4363 auroral lines. The differ- 
ent symbols correspond to the galaxies identified in the legend. References 
for the data are given in the text. The four Hll regions in NGC 4258 with 
a measured [O III] A4363 line are shown by the orange squares. The curves 
represent the upper and lower branch of the R23 indicator, obtained from a 
polynomial fit to the data. 

range and sufficiently strong to allow a good measurement 
of its intensity. This line, in combination with the much 
stronger [N II] A6583, can be used to measure the electron 
temperature in the low-ionization zone of an Hll region, 
in order to determine, for ex ample, the N + and + ionic 
abundances (Bresolin et al. 2004}. Here we can use it to infer 
the strength of the [N II] A6583 line, making a reasonable 
assumption for the value of the electron temperature. Then, 
the N2 = log([Nll] A6583/Ha) index can be used in the 
traditional way to assign the objects to the lower or upper R23 
branch. For example, assuming T e = 10,000 K (T e = 10,600 K 
for object 34, as derived from the [O III] A4363 line), 
from the temden routine in IRAF we obtain N2 = -0.72± 
0.08,-0.61 ±0.05, -0.77 ±0.05, and -0.94 ±0.06 for the 
four objects, respectively. It is generally adopted that for 
N2>-1.2, the upper branch solution o f R23 should be 
preferred (e.g. Kewley & Ellison 2008 1. The boundary 
between upper and lower branch could be safely lowered to 
N2 ~— 1.4, if we consider that the upper and lower branches 
intersect around 12 + log(O/H) = 8.0, when the abundances 
are determined empirically (see Fig. [3J. A plot of empirical 
([Olll] A4363-based) abundance s vs. N2 shows tha t this 
value corresponds to N2 = -1.5 (Pettini & Pagel 2004}. In 
any case, the N2 values we estimated above place the targets 
safely in the upper branch. Only once we approach a T e value 
of 14,000 K would the N2 estimates be below -1 .2. However, 
such high temperatures can be reasonably excluded by 
pointing out that the directly measured electron temperature 
values, approximately between 11,000 K and 12,000 K, 
are found for objects located at the largest galactocentric 
distances in NGC 4258, while the electron temp erature is ex - 
pected to decrease towards smaller radii ( Bresolin et al. 2009| l. 



In conclusion, there are strong arguments to favor the se- 
lection of upper branch calibrations of R23 for the whole H II 
region sample analyzed in this paper. 

4. THE OXYGEN ABUNDANCE GRADIENT IN NGC 4258 

The adoption of the R23 upper branch, justified in the pre- 
vious section, yields the oxygen abundances pr esented in 
colum ns (10) and (11) of Table |3] o btained wi th the McGaugh 
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E 
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3.30 


ZKH 5 


F 


12 18 55.470 


47 16 47.82 


7.27 


OK (-035-086) 


G 


12 18 59.308 


47 16 44.46 


4.47 


ZKH 7 


H 


12 18 57.856 


47 16 07.69 


7.84 


54C 
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12 18 58.614 
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8.00 


58C 
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12 18 58.997 


47 15 50.89 


8.06 


59C 
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K 


12 19 00.116 


47 15 36.95 


8.17 


69C 


L 


12 19 01.438 


47 15 25.66 


8.07 


OK (+024-170), 74C 



(1991] =M91) andl Pilyugin & thuan| ( pOOS] =P05) calibra- 
tions, respectively. For the derivation or the abundance gradi- 



Note. — Units of right ascension are hours, minutes and seconds, and units of dec- 
lination are degrees, arcminutes and arcseconds. Col. (1): Hll region identification. 
Col. (2): Right Ascension. Col. (3): Declination. Col. (4): Galactocentric distance 
in kpc . Col. (S ): Identification from Oey & Kennicutt 1 1993 =OK), Zarit sky et"aL] 
|1994| =ZK H), |Bresoiin et al.|(l999| =B) and |Courtes et al.||l993| studied by |Diaz| 
|et al.|2000) . 



ent the line fluxes of 12 additional H H regio ns publis hed by 
Oey & Kennicutt| ( 11993| l, |Zaritsky et al.| ( |1994| l and |Diaz et aL] 
(2000) were added to our sample. These objects are identified 
with letters in Fig. [T] and their positions and galactocentric 
distances can be found in Table [5] When the same target was 
observed by these authors, the average [O II] and [O III] line 
fluxes were computed. For objects from our sample that are 
in common with these authors (see Tableland Fig.|2|, only 
our new measurements were considered. 

The galactocentric O/H gradient in NGC 4258 is presented 
in Fig. [4] In this plot abundances for the Gemini H II region 
sample calculated with the M91 and P05 R23 calibrations are 
shown as full squares and full triangles, respectively. Open 
symbols (squares and triangles for M91 and P05, respectively) 
are used for the H II regions extracted from the literature. In 
the case of the abundances based on the P05 method objects 
with values of the excitation P < 0.3 were removed, since 
the calibration provided by P05 only covers P > 0.3, and we 
found that for objects with smaller excitation values the abun- 
dances could be overestimated. 

For the calculation of the radial abundance gradient a least- 
squares fit that accounts for errors in both coordinates (abun- 
dance and radius) was carried out. Because of the high incli- 
nation of the disk of NGC 4258 relative to the line of sight 
(78°), even a relatively small uncertainty in this angle can re- 
sult into an appreciable error in the computed deprojected dis- 
tances, especially for objects located near the minor axis. The 
source of the geometric parameters of the disk of NGC 4258 
adopted here, |van Albada (1980}, does not explicitly provide 
uncertainty estimates, but implies that the errors in the radio- 
determined position angle of the major axis and inclination 
are much smaller than 2° and 1°, respectively. A 1° uncer- 
tainty was assumed for both angles. We also accounted for 
an astrometric uncertainty of 0"4 in the H II region positions. 
The resulting error in the galactocentric distance reaches a 
maximum value of 6%, but for the majority of the targets it 
is on the order of 1-2%. 
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Figure 4. Galactocentric oxygen abundance gradient in NGC 4258. The abundances for the Gemini Hn region sample calculated with the M91 and P05 
calibrations are shown as full squares and full triangles, respectively. Open symbols (squares and triangles for M91 and P05, respectively) are used for the Hll 
regions extracted from the literature. The corresponding regression lines, calculated accounting for errors in both coordinates, are shown by the grey and green 
lines. The open circles show the oxygen abundance of the four H II regions where the [O III] A4363 was detected. In several cases the errors in the deprojected 
galactocentric distances are smaller then the size of the symbols used. 



With the oxygen abundances derived from the M91 calibra- 
tion the regression to the data (grey line in Fig.|4]i yields: 

12 + log(0/H) = 8.87 (±0.02) - 0.012 (±0.002) R kpc (1) 

where is the galactocentric distance measured in kpc. 
With the adopted values for the distance (7.2 Mpc) and 
the isophotal radius R25 (9'.31), the gradient slope can also 
be expressed as -0.23 ± 0.04 dexfl 25 -1 or -0.025 ± 0.004 
dexarcmin" 1 . 

In the case of the abundances derived from the P method of 
P05 the linear fit is (green line in Fig. [4]): 

12 + log(0/H) = 8.49 (±0.02) - 0.011 (±0.002) R kpc (2) 

corresponding to -0.21 ±0.03 dexfl 25 -1 or -0.023 ± 0.003 
dexarcmin" 1 . The dotted lines in Fig. show the 1 a scatter 
around these fits (0.08 and 0.07 dex, respectively). 

It is worth pointing out the following: 

- the slopes of the two fits agree within the uncertainties, but 
the intercepts differ by 0.38 dex. This is due to the afore- 
mentioned systematic offsets between different calibrations 
of strong-line indicators, in particular between one calibrated 
from photoionization models (M91) and one calibrated 
empirically (P05). 

- a look at Fig. [4] suggests that in the inner 8 kpc the gradient 
is virtually flat, with constant 12 + log(0/H) = 8.79 ± 0.07 
(M91) and 12 + log(0/H) = 8.40 ± 0.07 (P05), however this 
solution is not found to be statistically preferable over the 
single slope solution. 

- the result on the abundance gradient slope appears to be 
relatively robust. We considered the R23 calibrations by 



Zarits kyeTaL] (J1994j l and JTremonti et aD p004) , both of 
which are based on theoretical models, and obtained results 
that are in agreement with the two previous determinations, 
as Table |6| shows. It should also be noted that the Zaritsky 
|et al.| ( |1994| l calibration yields the largest intercept value, 
0.53 dex above the P05 calibration and 0.15 dex above the 
M91 calibration. 



The data points corresponding to the [O III] A4363-based 
abundances are shown as red circles (labeled A4363) in Fig. [4] 
Their positions in the plot are within ~ 1 a of the regression 
line for the abundances of the full sample based on the P05 
method. This relatively good agreement is not surprising, 
since the P method was calibrated empirically from H II re- 
gions whose abundances were obtained from the use of the 
auroral lines. For Cepheid-related work it is worth pointing 
out that the abundances tied to the [O III] A4363 detections 
(i.e. those obtained from the P05 calibr ation) are to be con- 
sidered in the r e -based scale (Ke nnicutt et al. 2003 |Sakai| 



eTaLl|2 004). Although the overall metallicity of NGC 4258 
in the T e scale might seem to be low, its characteristic oxy- 
gen abundance [12 + log(0/H) = 8.40 at R = 0AR 2 5] is nor- 
mal for this galaxy's luminosity and rotational velocity , ac- 
cording to the trends found by Pilyugin et al. ( |2004| l us- 
ing abundances of galaxies obtained from the P method, al- 
beit with an O/H ratio slightly below average. We can also 
compare NGC 4258 to a galaxy with a well-determined, T e - 
based abundance gradient like M101, which has an abso- 
lute luminosity M B = -20.9, similar to that of NGC 4258 
(M B = -20.8). The slope of the abundance gradient in M 101 
is -0.032 dex kpc" 1 , with a n extrapolated c entral abundance 
12 + log(Q/H) = 8.75 ± 0.05 ( |Bresolin|2007l ). The - 0.26 dex 



lower intercept for NGC 4258 can be understood by recalling 
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Table 6 

Results of linear regression to the abundance gradient 



R23 calibration 


Intercept 




Slope 








dexkpc 1 


dex/? 25 -1 


dex arcmin ' 


(1) 


(2) 


(3) 


(4) 


(5) 



McGaugh(1991l 



2005 1 



This paper 

8.87 ±0.02 -0.012 ±0.002 

8.49 ±0.02 -0.011 ±0.002 

9.02 ±0.03 -0.013 ±0.002 

8.94 ±0.02 -0.011 ±0.002 

Published 



Authors 

Zaritsky et"ai]|1994) 
Dutil & Koy|jr999^ 
Pilyugin eTffll^iUUq) 
Bono et al.ftZOOSf 



9.17 ±0.06 
8.86 ±0.02 

8.57 

8.55 



-0.029 ± 0.005 
-0.012 ±0.002 

-0.010 

-0.010 



-0.23 ±0.04 
-0.21 ±0.03 
-0.25 ±0.05 
-0.22 ±0.04 



-0.58 ±0.09 
-0.23 ±0.04 

-0.20 

-0.20 



-0.025 ± 0.004 
-0.023 ± 0.003 
-0.027 ± 0.005 
-0.023 ± 0.005 



-0.062 ±0.010 
-0.025 ± 0.004 

-0.021 

-0.021 



that strongly barred galaxie^] in addition to shallower gra- 
dients, also display significantly lower (0.3-0.4 dex) central 
abund ances compared t o non-barred galaxies (|Dutil & Roy| 
1999] see also Fig. 11 of |Vila-Costas & Edmunds) 1992) . 

As already pointed out, the absolute values of the chemical 
abundances are greatly dependent on the choice of abundance 
determination method. The P method by P05 provides results 
which are in good agreement with the auroral line-based 
abundances determined for a handful of objects. However, 
comparisons with independent methods would be highly 
desirable, in order to establish the absolute abundance 
scale in NGC 4258. This is relevant, for example, for the 
amount of the metallicity correction to the Cepheid-based 
distance modulus. Future spectroscopic studies of the stellar 
abundances in this galaxy are already planned by the author 
and his collaborators with observations at the Keck tel escope. 
It is also int eresting to note that the recent work by Konami 
et al. (2009) obtained a metallicity of the interstellar medium 
in NGC 4258 from observations wit h the Suzaku X-ray 
satellit e. Their result, ^0.5 Z in the |Anders & Grevesse" 
(1989) solar metallicity scale, corresponds approximately to 
12 + log(0/H) = 8.6, i.e. an intermediate value between those 
indicated by the two methods used here. 



5. PREVIOUS DETERMINATIONS OF THE ABUNDANCE 
GRADIENT 

The oxygen abundance gradient of NGC 4258 was already 
derived as part of earlier spectroscopic studies of larg e sam- 
ples of H II regio ns in different galaxies, in particular by |Zarit-| 
sky et al.| (fr994 = zk h), who also inc l uded s pectrophotomet- 
ric data obtained by |Oey & Kennicutt| ( |1993[l , for a total num- 
ber of about 15 objects. However, before a direct comparison 
with our results can be made it is necessary to homogenize 
the distances and diameters adopted by the different authors. 
In particul ar, ZKH adopted an isophotal radius of 7'92 from 
the RC2 ( |de Vau couleurs ~et al.||1976| ), which differs signifi- 
cantly from the value of 913 1 that we adopted from the RC3 
( de Vaucouleurs et al.|1991 



Thus, the gradient published by 



ZKH: 



12 + log(0/H) = 8.97 (±0.06) - 0.49 (±0.08) (R/R 25 -0.4) 

3 Evidence for the presence of a bar in NGC 4258 has been provided by 
van Albada 1 1980V 



becomes, after renormalization: 



12 + log(0/H) = 9.166 (±0.06) - 0.029 (±0.005) R kpc 

This is a significantly steeper gradient than the one we have 
derived earlier in Eq. ([T} and It appears that the ZKH re- 
sult is highly dependent on the abundance calculated for one 
of the two outermost H II regions, number 8 in their list (our 
number 30, for which we also obtained a [O III] A4363 detec- 
tion). The problem is not in the line fluxes, since we measure 
very similar [O II] A3727 and [O III] AA4959, 5007 line inten- 
sities. However, the application of ZKH's own calibration of 
R23 provides quite a low oxygen abundance for this partic- 
ular Hll region, 12 + log(0/H) = 8.41, about 0.4 dex lower 
than other objects at a similar galactocentric distance. This 
anomaly appears to be related to the high R23 value. In fact, 
two additional H II regions in our sample (24 and 36), which 
also have log(R23) > 0.9, display pecu liarly low O/H abun- 
dances when adopting the ZKH and the Tre monti et al.| ( |2004] > 
R23 calibrations. This suggests a systematic problem with 
the functional form of these calibrations at high R23 values. 
Therefore, we have excluded these three H II regions to calcu- 
late the regressions whose parameters are shown in Table [6] 

As part of an abundance study of eight early-type spi- 
ral galaxies, Dutil & Roy ( 1999| > measured the abundance 
gradient of NGC 4258, using a spectrophotometric tech- 
nique based on narrow-band imaging through filters cen- 
tered on important nebular emission lines and applied to 122 
Hll regions. They adopted as their abundance indicator the 
[O III] A5007/[N II] A6583 line ratio, calibrated as a function 



of O/H by |Edmunds & Page! (T1984). These authors adopted a 
distance to NGC 4258 of 7.3 Mpc, i.e. virtually the same one 
used here. They obtained the following oxygen abundance 
gradient: 



12 + log(0/H) = 8.86 (±0.02) - 0.012 (±0.002) R kpc 

which is the sa me we derive f rom our sample using the M91 
R23 calibration. |Dutil & Roy| ( |1999| l do not provide a list of 
line fluxes for their targets, so that we cannot compare their 
measurements with ours on an object-by-object basis. How- 
ever, the excellent agreement between the two gradient de- 
terminations i ndicates that the imaging spectrophotometry of 
|Dutil & Roy| ( |1999| l is not affected by important systematic 
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uncertainties, which could originate especially from the sub- 
traction of the stellar underlying continuum, and which could 
introduce b iases in the che mical abundances derived with this 
technique (|Dutil & Roy|2001|l. 



More recently, B ono et al.| (EDDSjl reconsidered the abun- 
dance gradient in NGC 4258 in their reanalysi s of the Cepheid 
variables discovered by |Macri et al.| ( [2006| . |Bono et aT7| ( |2008) l 
found that the comparison of the data with pulsation mod- 
els suggested both a lower metallicity and a shallower abun- 
dance gradient in NGC 4258 than one obtains from t he ZKH 
paper. Taking the few data points published by Dia z et aT] 
(2000), these authors derived an abundance gradient which, 
after renormalization to the quantities adopted in our paper, 
has a slope of -0.010 dexkpc" 1 , with a central abundance 
12 + log(Q/H ) = 8.55 . This agrees also with the study by |Pi-| 
|lyugin et al. (2004i, who used the P parameter method to 
measure abundance gradients fr om data in the lit erature for 
more than 50 galaxies. Finally, [Riess et aL] ( |2009a| l obtained 
abundances from Keck observations. They do not provide in- 
formation about their targets and their derived abundances, 
but from their Table 12 we infer a slope of -0.017 dexkpc" 1 . 
This gradient is possibly too steep because of the presence of 
two H II regions at large galactocentric distance and low O/H 
ratio (as judged from their Fig. 11), which might suffer from 
the same large-R23 issue mentioned earlier. 

To summarize, the evidence from most of the results pub- 
lished in the literature and from the analysis of the new 
Gemini data presented here is that the abundance gradient 
in NGC 4258 is very shallow, with a slope of approximately 
-0.012 dexkpc" 1 , in contrast with the steeper (by a factor of 
2.5) slope obtained by ZKH. This shallow abundance gradi- 
ent is not surprising, in view of the fact that barred galax- 
ies (like NGC 4258) and spiral galaxies of early type have 
smaller slopes compared to non-barred and late-type galax - 
ies (e.g. |Vila-C6stas" & Edmun dsp99^1 [Zaritsky et al.|1994) . 
For example, |Vila-Costas & Edmunds| ( |1992[ i found mat the 
barred galaxies in their sample have slopes that are shallower 
than -0.05 dexkpc" 1 , while non-barred galaxies can have 
slopes up to four times steeper. This behavior is als o pre- 
dicted by n umerical simulations of barred galaxies (Friedli 
|et al.|1994) . 

6. THE METALLICITY DEPENDENCE OF THE CEPHEID 
PERIOD-LUMINOSITY RELATION IN NGC 4258 



[Macri et al.| (|2006 ) have identified 281 Cepheid variables in 
two HST ACS fields of NGC 4258, whose centers are located 
179" (6.25 kpc, 'inner' field) and 490" (17.11 kpc, 'outer' 
field) from the galaxy nucleus (deprojected distances). Our 
spectroscopic coverage of H II regions in NGC 4258 was de- 
signed to allow the measurement of nebular abundances at 
the same galactocentric distances of the Cepheids, as visu- 
alized in Fig. [5] The extinction-corrected distance mo duli of 
the two fiel ds relative to the LMC, as determined by |Macri 
et al. (2006) from the BVI P-L relations of a restricted sample 
of Cepheids, are A/i„ = 10.87 (inner field) and A/x = 10.71. 
These authors ascribed this result to the difference in metal- 
licity between the two fields, as a consequence of the de- 
pendence of the P-L relation on metal content. Adopting 
the radial abundance gradient measured by ZKH, they fitted 
the distribution of the distance moduli derived for individ- 
ual Cepheids as a function of the O/H ratio, to obtain a value 
of the parameter 7 = 6fio/S\ogZ = -0.29 ± 0.09 (random) ± 
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Figure 5. The distribution of Cepheids in NGC 4258 I Macri et al. 2006] blue 
crosses) compared with the location of the H II regions studied in this work 
(circles for the new Gemini sample, squares for the objects drawn from the 
literature). The red circles represent the H II regions with the [O III] A4363 
de tections. The big c rosses show the centers of the two ACS fields observed 
by |Macri et al.|(2"006) . S-band image courtesy of L. van Zee. 



0.05 (systematic) magdex" 1 . This 7 value is in good agree- 
ment with other empirical determinations (see a re view of the 
exten sive data present in the literature by Romaniello et al. 
2008 ). Here, we just mention some of the recent results based 



on the differential analysis of Cepheids loc ated at different 
distanc es from the center of their host galaxy. Kennicutt et al. 
(1998) determined a value of 7 = -0.24 ±0.16 magdex" 1 
from HST observations of two fields in M101. More recently, 
|Scowcroft et al.| ( |2009| l performed a similar experiment in the 
galaxy M33, obtaining 7 = -0.29 ± 0. 1 1 mag dex" 1 . It is well- 
known that these empirical results, showing that metal-rich 
Cepheids are brighter than metal-poor ones, are at variance 
with the predictions from pulsation models, i.e. that Cepheids 
become fainter with increasing metallici ty, with a dependenc e 
on the passbands of the observations ( |Fiorentino et al.| 2007). 
The latter finding is backed by recent investigations o f the 
iron content of Galactic and Magellanic Cloud Cepheids ( |Ro-| 
Imaniello et al.|2008| ). f\ 

For the case of NGC 4258, as demonst rated in Sect . |4| and 
[5] the ZKH abundance gradient adopted by Macri et al. ( 2006 ) 
is much steeper than our new observations and other data in 
the literature indicate. Instead of a difference in oxygen abun- 
dance of <51og(0/H) = 0.31 dex between the inner and outer 
fields, as one would derive from the ZKH gradient, we mea- 
sure <51og(0/H) ~ 0.13 dex. Correspondingly the metallicity 
correction factor would become 7 = -0.69 magdex" 1 , an un- 

4 The 7 values reported here refer to distances measured from VI pho- 
tometry. The metallicity effect is predicted to vary with the use of different 
passbands. 
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Figure 6. Log(R.23) vs. 12 + log(0/H) for the same sample of extragalactic 
H II regions sh own in Fig. [3] T he curves in the upper part represent the R23 
calibrations by Zaritsky ef~al. 1 1994 continuous line, ZKH) and Tremonti 
|etlu~l(2004] dashed line, T04). 

acceptably large value from both the empirical and theo retical 
standpoints. We are drawn to conclude, together with |Bono| 
[etaLl ( [2008) l, that, because of its shallow gradient, NGC 4258 
(like any other barred galaxy) is not the ideal laboratory to 
perform differential studies of the metallicity dependence of 
the Cepheid pulsation properties. 

What is then the origin of the observed distance modulus 
difference between the inner and outer fields in NGC 4258? 
While we cannot presently discard the effects of varying 
Cep heid helium abundance or differences in the reddening 
law ( |Mager et al.|2008] l, the presence of systematic errors af- 
fecting the Cephe id distance of the o uter field appea rs as a 
possible solu tion ( |di Ben edetto 2008). In particular, Mager 
|et al.| ( |2008| l showed how the outer field Cepheids, contrary 
to the inner field ones, do not evenly populate the instability 
strip, projected into the P-L relation. Instead, their magni- 
tudes are skewed towards faint levels, thus artificially increas- 
ing the distance modulus. 

6.1. The metallicity of NGC 4258 on the T e scale 

The detection of the [O III] A4363 line and the derivation of 
O/H ratios from the strength of this line in four of the outer 
H II regions has allowed us to place the chemical abundance 
of NGC 4258 on the T e scale on rather solid evidence, in com- 
bination with the use of the P05 calibration of the R23 metal- 
licity indicator for the innermost part of the galaxy. From 
Eq.|2]we obtain: 

inner field: 12 + log(0/H) = 8.42 ± 0.03 
outer field: 12 + log(0/H) = 8.30 ±0.06 

If for the outer field we take the A4363-based abundance for 
objects 24 and 34 (their galactocentric distance matches that 
of the center of the outer ACS field) we obtain 12 + log(0/H) 
= 8.23 ±0.09. Given the size of the errors, and the statisti- 
cal nature of the R23 indicator, we do not regard this small 
difference as a significant discrepancy, and adopt the former 
result for consistency with the inner field abundance, which 
can only be derived from the linear regression. 

It is interesting to compare these abundances with those 
proposed by other authors working on the Cepheid distance 
scale. Sak ai et aL"1 (|2"004) obtained a transformation between 
nebular oxygen abundances of Cepheid host galaxies based 
on the ZKH abundance scale, and those one derives using the 
auroral line method, which is based on the ability to directly 
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Figure 7. Comparison between the oxygen abundances in the ZKH and T e 
(i.e. [O III] A4363-based) scales. The sample of H II regions is the same as 
in Fig. [3] The linear regression to the H II region sample is shown by the 
straightTine, and the ±1<j scatter is represented by the shaded area. The 
high-metallicity (R23 upper branch) galaxy data points of Sakai et al. ( 2004 ) 
ar e shown by the blu e triangles. The dashed line is the linear ht to the latter 
by|Saha et alT]j2006). 



measure the electron temperature of the ionized gas (hence 
the T e -based' definition). As already stressed, the T e abun- 
dances are systematically lower than those derived from the 
ZKH and most of the other strong-line metallicity indicators 
( |Kennicutt et~aLl[2003l |Bresolin et ail|2004"l ). The effect is 
stronger at high metallicity, thus compressing the metallicity 
range spanned by star-forming galaxies. This can be appreci- 
ated in Fig. [6] where the same data points of Fig [3] are shown, 
together with th e curves corres ponding to the ZKH (contin- 
uous line) and Tre monti et al.| ( |2004| dashed line) R23 cali- 
brations. These two curves are seen to be diverging from the 
r c -based abundances, shown by the dots, as the metallicity in- 
creases, with the effect being stronger in this example for the 
ZKH calibration. 

It should also be recalled that most of the nebular work 
of the past few years has adopted the recent solar oxy- 
gen abundance derived from 3-D hydrodynamic al models of 
the Sun, which yield 12 + log(0/H ) Q = 8.69 ( |Allende Prieto 
eTlLp OOl 1 Asplun d~et~aT1[2U09] l, 0.24 dex below the |A~ 
ders & Grevesse| ( |1989| l value. In this scale, for example, 
the LMC is found to have a nebular oxygen abundance of 
12 + log(O/H) LM c = 8.36±0.10. This value is the mean ob- 



tained by the author (see Appendix) from a reanalysis of 
published emission line fluxes carried out with consistent 
atomic data. This is in excellent agreement with results from 
other stellar indicator s, for example Cepheids, for which Ro- 
maniel FetaLl p008 ) found [Fe/H] LMC =-0.33±0.13, corre- 
sponding to 12 + log(O/H) = 8.36±0.13. 



Saha et al. (2006) and Bono et al. 



( 2010| l performed poly- 



akai et alTf (2004) data points, producing 
O/H abundance ratios in the ZKH system 



nomial fits to the 

a relation between O/H abundance ratios in the ZKH system 
and th e T e system. Using the data in Table 7 of Bo no et aL] 
( |2010| > we find that in the case of NGC 4258 the transformed 
12 + log(0/H) abundances are 0.2 dex larger than the abun- 
dances that we measure: 8.64 vs. 8.42 for the inner field, 8.50 
vs. 8.30 for the outer field. To investigate the origin of this dis- 
crepancy, the relation between the oxygen abundances in the 
ZKH and T e scales for the same sample of H II regions shown 
in Fig. [3] was considered. The result is i llustrated in Fig. [7 
where the abundances from Table 4 of |Sakai et al.| (|2004 



open triangles) and the linear fit to the latter by Saha et al. 



(2006, dashed line) are included. Only points in the R23 upper 
branch are shown, because the ZKH calibration only applies to 
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that regime (the transformation between the two abundance 
scales below 12 + log(0/H)zKH — 8.2 is thus undefined). The 
linear fit to the individual H II region data (con tinuous line) in 
Fig. |T|lies 0.15 dex below the Saha et al. (2006) regression at 
x = 12 + log(0/H) ZKH = 8.5. The new transformation is: 



12 + log(0/H) Te = 0.69 (±0.06) x + 2.30 (±0.51). (3) 



Only two data points from Sakai et al. (2004), the inner and 
outer fields of M101, fall within the l a scatter (±0.11 dex) 
of this regr ession. The orig in of the r e -based abundances pro- 
vided by Sakai et al. (2004, Table 4) is not specified in their 
paper, so it is difficult to ascertain the reason for the discrep- 
ancy for the remaining upper branch metallicities. It should 
also be added that at the time of their publication the system- 
atic offset between strong-line abundances and direct abun- 
dance s in Hll regions had just become clear (Kennicutt et al. 
2003), and that since then the number of high-metallicity H II 
regions with auroral line detections has increased significantly 
(e.g. |Bresolin et al.|2009) >. 

The comparison between our new and the previous ZKM-T e 
abundance transformations explains the reason why the metal- 
licities we derive for the two fields in NGC 42 58 are smaller 
than th ose obtained from the Saha et al. (2006) orBo noet al.| 
( |2010[ ) transformations. It also indicates that the latter should 
be revisited according to Eq. [3] if one wishes to transform in 
a consistent way the strong-line abundances in the ZKH scale 
to the T e scale. This is done in the Appendix for a number of 
Cepheid host galaxies. 

We would also like to point out that the fact that the mul- 
tiplicative coefficient forx in Eq. [3] differs significantly from 
unity (the same information is conveyed by the divergence of 
the ZKH curve in Fig. [6] from the 7^-based empirical abun- 
dances with increasing metallicity) implies that abundance 
gradient slopes can depend on the choice of abundance de- 
termination methods. The effect is small in the case of 
NGC 4258, where the gradient slope is very shallow and the 
metallicity range over the full extent of the galactic disk is 
~ 0.2 dex. However, for galaxies with steeper abundance gra- 
dients the effect could become significant. 

6.2. Consequences for the Cepheid distances 

The abundances in the T e scale place the inner field of 
NGC 4258 essentially at the metallicity of the LMC, rather 
than at the Galactic value (12 + log(0/H)~8.7), which has 
been usually adopted in the lit erature to estimate the distance 
modulus to NGC 4258 (e.g. jBenedict et al.Jj2007] ). How 
would this affect the C epheid! distance to NGC 4258 from 
the |Macri et~aL ( 2006 1 data? To answer this question, two 
methods have been considered: 

a) we take the extinction-corrected distance modulus of the 
inner fi eld of NGC 4258 relative to the LMC from lMacri et al.l 
(2006), A^„ = 10.71 ±0.04 r ± 0.05. s , and apply no metallicity 
correction, since the O/H abundance is the sa me. Thus, 
adopti ng /zl mc = 18.5 from the HST Key Project (Freedman 
[etaLl [2001) , we obtain a*ngc 4258 =29.21 ± 0.o4 r ± 0.05 s . 
Alternatively, if we assume the maser distance modulus 
to NGC 4258, 29.29 ± 0.09 r ± 0.12 s we would obtain 
Mlmc= 18.58±0.10,.±0.13 s . 



days, for consistency with the procedure adopted by the Key 
Project. For the inner field they obtain /io = 29.26 ±0.03. 
The metallicity correction applied by the Key Project is 
8[i z = -0.2[(O/H) - (0/H)lmc] magdex" 1 . Since we are 
testing the idea that the inner field has the same metallicity 
of the LMC, this correction would be zero in our case. Thus, 
MOJVGC 4248 = 29.26 ± 0.03. 

We also note that the adoption of an LM C-like metallic- 
ity for the inner field explains the finding by Tammann et al. 
( 2008 ) that the BVI P-L relations in this field agree with those 
in the LMC, rather than those for the higher Galactic metal- 
licity, without the need to invoke the existence of a second 
parameter, such as the helium abundance. In virtue of the 
shallow abundance gradient, the outer field is also at a simi- 
lar metallicity. We can also conclude that the adoption of the 
T e abundance scale for NGC 4258, as presented here, would 
yield results that are compatible with the Key Project and with 
the geometric distance from the water masers. 

6.3. Reliability of the absolute T e abundances 

Although our analysis of H II region data can be used to 
place the metallicity of the interstellar medium of NGC 4258 
on the T e scale, there are some important considerations that 
need to be taken into account when dealing with nebular abun- 
dances obtained from the [O III] A4363 auroral line. First is 
the effect of dust depletion onto dust grains. Recent esti- 
mates of the oxygen depletion factor obtained from the study 



of the Orion neb ula are in the range b etween 0.08 de x ( Es- 
teban et al.||2004] > and 0.12 dex ( jMesa-Delgado et al|2009| ). 
Secondly, accounting for the presence of temperature fluctua- 



Secondly, accounting for the presence of temperature 
tions in the ioni zed gas has the effect of increasing the metal 
abundances (Peimbert 1967), yielding results that are com- 
patible with those from metal recombination lines. For these 
reasons, the r e -based abundances are often regarded as lower 
limits to the real abundances. 

In the case of the Orion nebula, whose T e -bas ed oxy- 
gen abundan ce is 12 + log(0/H) = 8.51 ±0.03 ( |Esteban| 
|et al. | [20041 ), the sum of these two effects is about 0.25 
dex. Recently Simon-Diaz & Stasiriska (2010 1 have shown 
that the oxygen abundance for the gas phase (from metal 
recombination lines) in the Orion nebula, increased by the 
amount locked in dust grains, is consistent with the value of 
12 + log(0/H) = 8.74±0.04 derived from 13 B stars in the Ori 
OBI association. This compares very well with the Galac- 
tic neigh borhood value of 12 + l og(Q/H) = 8.76 ± 0.03, calcu- 
lated by Przybilla et al. (2008) from B stars. These results 
suggest that the r c -based results could systematically under- 
estimate the oxygen abundance by 0.25 dex. 

On the other hand, the addition of the depletion correction 
factor to the recombination line abundance of 30 Dor in the 



LMC gives an abundan ce, 12 + log(0/H) = 8.61 ±0.03 (Peim- 



bert & Peimbert 2010 ), which is clearly at variance with re- 
spect to the B-star [12 + log(Q/H) = 8.33 ± 0.08,|Hunter et al. 



Pi 
les 



2007 1 and Cepheid (Romaniello et al.|2008 ) metallicities, un- 



ess it is assumed that the present-day chemical composition 
of the LMC is highly non-uniform. Therefore, in this case 
the r e -based oxygen abundance is in excellent agreement with 



the stellar one, as is also found for the SMC (e.g. Trundle 



& LennonpoB) , where 12 + log(0/H) = 8.03 ± 0.09 (see Ap- 
pendix). 



b) we follow Mager et al. (2008 1, who restricted the Macri 



et al.| ( |2006| ) sample to Cepheids with periods above 10 
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6.4. The case ofM33 
We also br iefly consider the case of M33, where |Scowcroft| 



et al. 



(2009) have determined 7 = -0.29 ±0.11 mag dex" 
from the difference in distance modulus between an inner and 
an outer field, and adopting a metallicity gradient obtained 
from Hll regions which is very steep (-0.19 dex kpc" 1 ) in 
the central 3 kpc, and rathe r shallow (-0.038 dexkpc" 1 ) be- 
yond this radius (taken from Magr ini et al.|2007) . The choice 
of this bimodal abundance gradient is questionable, and the 
most rece nt works on the chemical composition of the ion- 
ized gas (|Crockett et al.|[2006l [Rosblowsky & SimonJI2008j 
|Bresolin et al.||2010| l and young supergiants ( |U et al.j|2009f 
in M33 yield a single-slope gradient, with a shallow value 
throughout the disk of the galaxy. The most recent and ar- 
guably one of the most accurate determinations to date of 
the H II region radial oxygen ab undance gradient in M33 has 
been obtained by Bresoli n et al.| ( |2010] l, who measured a slope 
of -0.030 ± 0.008 dexkpc -1 from [O III] A4363 detections, 
which extend into the inner kpc of the galaxy. They found 
no indication of an upward turn in the nebular abundances 
in the inner 3 kpc, where they measured 12 + log(Q/H) ~ 8.4, 



instead of 12 + log(0/H) = 8.85, as used by |Scowcroft et aL 
(2009). The adoption of such a shallow abundance slope 
yields a total range in O/H between the inner and outer regions 
of M33 of only 0.14dex. This, combined with the distance 
moduli obtained at the two different galactocentric distances 



by Scowcroft et al. (12009), would result in an exceedingly 
strong metallicity effect on the Cepheid-derived distances, 
with 7 ~ -1 .2 mag dex" 1 . 

We note that a bias from small number statistics, similar 
to the one encountered i n NGC 4258, could also be affect- 
ing the M33 result by |Scowcroft et al. ( |2009| ). As can be 
seen from the P-L and Period- Wesenheit relations shown in 
their Fig. 5 and 6, the scatter for the outer field Cepheids at 
constant period is much smaller than that for the inner field, 
indicating that the width of the instability strip is not suffi- 
ciently sampled. This effect is exacerbated by the scarcity 
of Cepheids with periods > 10 days in their outer field (only 
two are found). We reach the same conclusion we drew for 
NGC 4258, that 7 cannot be reliably derived in M33 with 
present Cepheid samples. A more reliable approach to con- 
strain the metallicity effect on Cepheid-based extragalactic 
distances is represented by the comparison with the distances 
provided by the TRGB method as a function of the host galaxy 
metallicity JSakai et al |[2004{ |Tammann et~aT1 [2008 ; Bono 
|et al.|20T0) ), provided the latter is measu red c onsistently, for 
example in the T e scale, as explained in § 6.1' 



7. SUMMARY 

This paper combined new deep spectroscopic observations 
of 36 H II regions in NGC 4258 obtained with the Gemini 
telescope with data from the literature, to measure the oxygen 
abundance gradient for this galaxy. Strong-line abundances 
have been derived fr om two different calibrations of t he R23 
indicator, the one by McGaugh|( |1991| ) and the one by |Pilyu-] 
|gin & Thuan ( |2005[ ). For four of the outermost H II regions 
the |OlllJA4363 auroral line has been detected, and used to 
measure the nebular abundances directly in the T e scale. The 
main results are summarized here: 

1. The radial abundance gradient in NGC 4258 has a slope 
of approximately -0.012 dexkpc" 1 , which is a fa ctor of 2.5 
shallower than the value reported by|Zantsky^et al. ( 1994} and 



adopted by Macri et al. ( 2006 ) to infer a metallicity correction 
of Cepheid distances of 7 = -0.29 mag dex" 1 . 

2. Adopting the shallow gradient resulting from the new anal- 
ysis would yield an unrealistically large metallicity correction 
for the Cepheid distances. This indicates that NGC 4258, with 
its present sets of Cepheid observations, is not suited for a dif- 
ferential metallicity study to measure 7 reliably. 

3. The [Olll]A4363 detections and the P05 R 23 calibration 
provide the oxygen abundances in the T e scale. In this scale, 
the inner field of NGC 4258 has an LMC-like O/H ratio, and 
not a Galactic one. 

4. The same systematic biases from small number statistics 
that might explain the difference in distan ce moduli between 
the inner and outer regions of NGC 4258 ( |Mager et al.|2008| ) 
se em to affect also the dif ferential analysis of M33 Cepheids 
by |Scowcroft et al.| ( [2009l ). 

5. Th e transformation betw een 'old' Hll region metallicities 
in the[Zaritsky et aT ( 1994 1 system and the T e system by Sakai 
|et al.| ( |2"0"04| ) is affected by a systematic error of about 0.2 dex. 
Using recent observations of H II regions with auroral line- 
based abundances, a revised transformation has been calcu- 
lated. 

The main results of this work can be tested in the near fu- 
ture, with the availability of larger samples of Cepheids in 
the outer regions of NGC 4258 and M33, to remove potential 
statistical biases, and with the derivation of blue supergiant 
metallicities, which will help to establish the absolute scale of 
the chemical abundances in NGC 4258. 



FB gratefully acknowledges the support from the National 
Science Foundation grants AST-0707911 and AST-1008798, 
and thanks R.P Kudritzki for discussions and comments on 
an earlier draft of this paper. The Gemini observations used 
in this paper were obtained under program GN-2010A-Q-25. 



Facility: GeminkGillett (GMOS) 



REFERENCES 



Allende Prieto, C, Lambert, D. L., & Asplund, M. 2001, ApJ, 556, L63 
An, D., Terndrup, D. M., & Pinsonneault, M. H. 2007, ApJ, 671, 1640 
Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 
Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 
Benedict, G. E, et al. 2007, AJ, 133, 1810 

Bono, G., Caputo, E, Fiorentino, G., Marconi, M., & Musella, I. 2008, ApJ, 
684, 102 

Bono, G., Caputo, E, Marconi, M., & Musella, I. 2010, ApJ, 715, 277 
Bresolin, F. 2007, ApJ, 656, 186 

Bresolin, F. 2008, in The Metal-Rich Universe, ed. G. Israelian & 
G. Meynet, 155 

Bresolin, E, Garnett, D. R., & Kennicutt, R. C. 2004, ApJ, 615, 228 
Bresolin, E, Gieren, W., Kudritzki, R., Pietrzyhski, G., Urbaneja, M. A., & 

Carraro, G. 2009, ApJ, 700, 309 
Bresolin, E, Kennicutt, Jr., R. C, & Garnett, D. R. 1999, ApJ, 510, 104 
Bresolin, E, Pietrzyriski, G., Urbaneja, M. A., Gieren, W., Kudritzki, R.-P, 

& Venn, K. A. 2006, ApJ, 648, 1007 
Bresolin, E, Schaerer, D., Gonzalez Delgado, R. M., & Stasihska, G. 2005, 

A&A,441,981 

Bresolin, E, Stasihska, G., Vi'lchez, J. M., Simon, J. D., & Rosolowsky, E. 

2010, MNRAS, 404, 1679 
Bresolin, E, Urbaneja, M. A., Gieren, W., Pietrzyhski, G., & Kudritzki, R.-P. 

2007, ApJ, 671,2028 



Abundances in NGC 4258 



13 



Courtes, G., Petit, H., Hua, C. X, Martin, P., Blecha, A., Huguenin, D., & 

Golay, M. 1993, A&A, 268, 419 
Crockett, N. R., Garnett, D. R., Massey, P., & Jacoby, G. 2006, ApJ, 637, 

741 

de Vaucouleurs, G, de Vaucouleurs, A., Corwin, Jr., H. G, Buta, R. J., 
Paturel, G., & Fouque, P. 1991, Third reference catalogue of bright 
galaxies (Springer- Verlag Berlin Heidelberg New York) 

de Vaucouleurs, G, de Vaucouleurs, A., & Corwin, J. R. 1976, Second 
reference catalogue of bright galaxies (University of Texas Press) 

di Benedetto, G. P. 2008, MNRAS, 390, 1762 

Diaz, A. I., Castellanos, M., Terlevich, E., & Luisa Garcia- Vargas, M. 2000, 

MNRAS, 318, 462 
Dufour, R. J. 1975, ApJ, 195, 315 
Dufour, R. J., & Harlow, W. V. 1977, ApJ, 216, 706 
Dufour, R. J., Shields, G. A., & Talbot, Jr., R. J. 1982, ApJ, 252, 461 
Dutil, Y., & Roy, J.-R. 1999, ApJ, 516, 62 
— . 2001, AJ, 122, 1644 

Edmunds, M. G, & Pagel, B. E. J. 1984, MNRAS, 211, 507 

Esteban, C, Bresolin, E, Peimbert, M., Garcfa-Rojas, J., Peimbert, A., & 

Mesa-Delgado, A. 2009, ApJ, 700, 654 
Esteban, C, Peimbert, M., Garcfa-Rojas, J., Ruiz, M. T, Peimbert, A., & 

Rodriguez, M. 2004, MNRAS, 355, 229 
Ferrarese, L., et al. 2000, ApJS, 128, 431 

Fiorentino, G, Marconi, M., Musella, I., & Caputo, F. 2007, A&A, 476, 863 
Freedman, W. L., & Madore, B. F. 2010, ARA&A, 48, 673 
Freedman, W. L., et al. 1994, ApJ, 427, 628 
— . 2001, ApJ, 553, 47 

Friedli, D., Benz, W., & Kennicutt, R. 1994, ApJ, 430, L105 
Garcia-Rojas, J., Esteban, C, Peimbert, A., Rodriguez, M., Peimbert, M., & 

Ruiz, M. T. 2007, Revista Mexicana de Astronomia y Astrofisica, 43, 3 
Garnett, D. R., Shields, G. A., Skillman, E. D., Sagan, S. P., & Dufour, R. J. 

1997, ApJ, 489, 63 
Guseva, N. G, Papaderos, P., Meyer, H. T., Izotov, Y. I., & Fricke, K. J. 

2009, A&A, 505, 63 
Herrnstein, J. R., Moran, J. M., Greenhill, L. J., & Trotter, A. S. 2005, ApJ, 

629,719 

Herrnstein, J. R., et al. 1999, Nature, 400, 539 

Hook, I. M., J0rgensen, I., Allington-Smith, J. R., Davies, R. L., Metcalfe, 
N., Murowinski, R. G., & Crampton, D. 2004, PASP, 116, 425 

Humphreys, E. M. L., Reid, M. J., Greenhill, L. J., Moran, J. M., & Argon, 
A. L. 2008, ApJ, 672, 800 

Hunter, I., et al. 2007, A&A, 466, 277 

Izotov, Y. I., Thuan, T. X., & Lipovetsky, V. A. 1994, ApJ, 435, 647 
Kennicutt, R. C, Bresolin, F, & Garnett, D. R. 2003, ApJ, 591, 801 
Kennicutt, Jr., R. C, Stetson, P. B., Saha, A., et al. 1998, ApJ, 498, 181 
Kewley, L. J., & Ellison, S. L. 2008, ApJ, 681, 1183 
Kniazev, A. Y, Grebel, E. K., Pustilnik, S. A., Pramskij, A. G, & Zucker, 

D. B. 2005, AJ, 130, 1558 
Kobulnicky, H. A., Skillman, E. D., Roy, J., Walsh, J. R., & Rosa, M. R. 

1997, ApJ, 477, 679 
Konami, S., et al. 2009, PASJ, 61, 941 

Kudritzki, R.-P, Urbaneja, M. A., Bresolin, F, Przybilla, N., Gieren, W., & 

Pietrzyhski, G. 2008, ApJ, 68 1 , 269 
Kuzio de Naray, R., McGaugh, S. S., & de Blok, W. J. G. 2004, MNRAS, 

355, 887 

Lee, H., Skillman, E. D., & Venn, K. A. 2005, ApJ, 620, 223 
— . 2006, ApJ, 642, 813 

Leonard, D. C, Kanbur, S. M., Ngeow, C. C, & Tanvir, N. R. 2003, ApJ, 
594, 247 



Lopez-Sanchez, A. R., Esteban, C, Garcfa-Rojas, J., Peimbert, M., & 

Rodriguez, M. 2007, ApJ, 656, 168 
Macri, L. M., Stanek, K. Z., Bersier, D., Greenhill, L. J., & Reid, M. J. 2006, 

ApJ, 652, 1133 

Mager, V. A., Madore, B. F, & Freedman, W. L. 2008, ApJ, 689, 721 
Magrini, L., Vflchez, J. M., Mampaso, A., Corradi, R. L. M., & Leisy, P. 

2007, A&A, 470, 865 
McCommas, L. P., Yoachim, P., Williams, B. F, Dalcanton, J. J., Davis, 

M. R., & Dolphin, A. E. 2009, AJ, 137, 4707 
McGaugh, S. S. 1991, ApJ, 380, 140 

Menzel, D. H, Aller, L. H, & Hebb, M. H. 1941, ApJ, 93, 230 
Mesa-Delgado, A., Esteban, C, Garcfa-Rojas, J., Luridiana, V, Bautista, M., 

Rodriguez, M., Lopez-Martfn, L., & Peimbert, M. 2009, MNRAS, 395, 

855 

Newman, J. A., Ferrarese, L., Stetson, P. B., Maoz, E., Zepf, S. E., Davis, 

M., Freedman, W. L., & Madore, B. F. 2001, ApJ, 553, 562 
Oey, M. S., & Kennicutt, Jr., R. C. 1993, ApJ, 41 1, 137 
Pagel, B. E. J., Edmunds, M. G, Blackwell, D. E., Chun, M. S., & Smith, G. 

1979, MNRAS, 189, 95 
Pagel, B. E. J., Edmunds, M. G, Fosbury, R. A. E., & Webster, B. L. 1978, 

MNRAS, 184, 569 
Pefia, M., Stasinska, G, & Richer, M. G. 2007, A&A, 476, 745 
Peimbert, A., & Peimbert, M. 2010, ApJ, 724, 791 
Peimbert, M. 1967, ApJ, 150, 825 
Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59 
Pilyugin, L. S., & Thuan, T. X. 2005, ApJ, 631, 231 
Pilyugin, L. S., Vflchez, J. M., & Contini, T. 2004, A&A, 425, 849 
Przybilla, N., Nieva, M.-F, & Butler, K. 2008, ApJ, 688, L103 
Riess, A. G, Li, W., Stetson, P. B., Filippenko, A. V, Jha, S., Kirshner, R. P., 

Challis, P. M., Garnavich, P. M., & Chornock, R. 2005, ApJ, 627, 579 
Riess, A. G, et al. 2009b, ApJ, 699, 539 
— . 2009a, ApJS, 183, 109 
Romaniello, M., et al. 2008, A&A, 488, 731 
Rosolowsky, E., & Simon, J. D. 2008, ApJ, 675, 1213 
Russell, S. C, & Dopita, M. A. 1990, ApJS, 74, 93 

Saha, A., Thim, R, Tammann, G. A., Reindl, B., & Sandage, A. 2006, ApJS, 
165, 108 

Sakai, S., Ferrarese, L., Kennicutt, R. C, & Saha, A. 2004, ApJ, 608, 42 
Scowcroft, V, Bersier, D., Mould, J. R., & Wood, P. R. 2009, MNRAS, 396, 
1287 

Seaton, M. J. 1979, MNRAS, 187, 73P 

Simdn-Dfaz, S., & Stasinska, G. 2010, ArXiv e-prints 

Stanghellini, L., Magrini, L., Villaver, E., & Galli, D. 2010, A&A, 521, A3 

Tammann, G. A., Sandage, A., & Reindl, B. 2008, ApJ, 679, 52 

Thim, R, Tammann, G. A., Saha, A., Dolphin, A., Sandage, A., Tolstoy, E., 

& Labhardt, L. 2003, ApJ, 590, 256 
Tremonti, C. A., et al. 2004, ApJ, 613, 898 

Trundle, C, Dufton, P. L., Lennon, D. J., Smartt, S. J., & Urbaneja, M. A. 

2002, A&A, 395, 519 
Trundle, C, & Lennon, D. J. 2005, A&A, 434, 677 

Tsamis, Y. G, Barlow, M. J., Liu, X., Danziger, I. J., & Storey, P. J. 2003, 

MNRAS, 338, 687 
Turner, D. G. 2010, Ap&SS, 326, 219 

U, V, Urbaneja, M. A., Kudritzki, R.-P., Jacobs, B. A., Bresolin, R, & 

Przybilla, N. 2009, ApJ, 704, 1 120 
van Albada, G. D. 1980, A&A, 90, 123 

van Leeuwen, R, Feast, M. W., Whitelock, P. A., & Laney, C. D. 2007, 

MNRAS, 379, 723 
Vila-Costas, M. B., & Edmunds, M. G. 1992, MNRAS, 259, 121 
Webster, B. L., & Smith, M. G. 1983, MNRAS, 204, 743 
Zaritsky, D., Kennicutt, Jr., R. C, & Huchra, J. P. 1994, ApJ, 420, 87 



APPENDIX 



A. METALLICITIES OF CEPHEID HOST GALAXIES IN THE T e SCALE 

In this Appendix the oxygen abundances of Cepheid host galaxies are calculated self-consistently, using (a) observati ons of 
the [O III] A4363 auroral line or (b) the conversion given in Eq. |3]between abundances obtained from the Zaritsky et al. ( 1994 



ZKH) calibration of R23 and r e -based abundances. This equation was obtained by comparing upper branch abundances for a 
large sample of H II regions measured using the two methods. It should be recalled that the ZKH calibration only holds for the 
metal-rich branch of R23, and therefore the transformation is not defined for the low-metallicity branch (12 + log(0/H)7> < 8). 
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Table 7 

Cepheid host irregular galaxies with [O III] A4363 detections 



Galaxy 

LMC 

SMC 

IC 1613... 

WLM 

NGC3109. 
NGC 6822 . 
Sextans A . 
Sextans B . . 
NGC 5253 . 



12 + log(Q/H) 



Sources of emission line data 



8.36 ±0.10 
8.03 ± 0.09 

7.78 ± 0.05 
7.82 ± 0.08 

7.79 ± 0.08 
8.12 ±0.07 
7.48 ± 0.02 
7.59 ± 0.20 
8.20 ± 0.03 



[Dufou r ( 1975 1; Pagel et al. 1 1978 1; Dufour 7t^(l982) 
|Kussell 8rT3op i^(llMi);| fe 
Dutour I 19/5 1; Dutour & Harlow 1 19/7 1; Pagel et al. ( 1978 1 
jkussell teT)opitaj|T990| ; Isamls eTaT^OOSJ 



[Bresolin et al.||2 U07 
|Lee et al.|pTO5 
Ifena et aLlUUU/j 



1 samis et ai 



|Lee et al.|f20Ubf 
|Kniazev et al.H2005 
|isjiiazev eTaT|j/u iD : 

|Kobulnicky~eTal.|{1997} ; Lopez-Sanchez et al. {2007} 



Table 8 

Cepheid host spiral galaxies with [O III] A4363 detections 



Galaxy 12 + log(0/H)j- ( , Sources of emission line data 



M33 (inner) 1 8.39 ± 0.04 

M33 (outer) 1 8.26 ± 0.07 

MlOl (inner) 2 8.64 ± 0.06 

MlOl (outer) 2 8.26 ± 0.08 

NGC 300 8.33±0.13 

NGC 2403 8.28±0.11 


Bresolin et al. 
Bresolin et al. 
Kennicutt et al 


2010 
2010 
{ME 








Bresolin 


2007 


Kennicutt et al 


W: 




Bresolin 




Bresolin et al. 
Uarnett et al.jj" 


Uouy 
y9/| 




SmwCToft Rt al ,l|?.00q| 
Kennicutt et al.|fl998) 



A.l. Irregular and dwarf galaxies with abundances determined from [O III] X4363 detections. 

Emission line data for nine irregular and dwarf galaxies with [O III] A4363 line detections were drawn from the literature, and 
the oxygen abundan ces for individual H I I regions were calculated using the same method (iraf's temden and ionic routines) 
and atomic data (see Bre solin et al.|2009| for details) used for our T^based abundance determinations in NGC 4258. The average 
abundances and standard deviations are reported in Table [7] together with the literature sources for the emission line data. 



A.2. Spiral galaxies with abundance gradients determined from [O III] X4363 detections. 

Four spiral galaxies with a well-determined oxygen abundance gradient based on the detection of [O III] A43 63 are pres ented 
in Table [8] Cepheids in two fields at different galactocentric distances have been observed in M33 ( |Scowcroft et al.|2009| ) and 
MlOl Kennicut t et al.| ( [T998[ l. The abundances in the inner and outer fields of M33 and MlOl were calculated from the following 
expressions for the abundance gradients: 



M33: 12 + log(0/H) = 8.42 (±0.03) - 0.26 (±0.07) R/R 25 , with R 25 = 35'A ( |Bresolin et al.|201Q| >. 

(Rinner = 3.6, Router = 21.1) 



Ml 01: 12 + log(0/H) = 8.75 (±0.05) - 0.90 (±0.07) R/R 25 , with R 25 = 14'A ( jBresolin^OOT) . 

(Rinner ~ 1 -7, Router — 7.9) 



The Cepheids studied in NGC 300 cover a wide range in radius, and one could crudely divide them into an inner and an outer 
sample, using the following expression for the abundance gradient: 



NGC 300: 12 + log(Q/H) = 8.57 (±0.02) - 0.41 (±0.03) R/R 25 , with R 25 = 9.'75 ( |Bresohn et 31^2009 i 



The oxygen abundance for NGC 300 in Table [8] is the average for H II regions at galactocentric distances large r than 0.3 R 2 s. 
where the bulk of the Cepheids is located. In the case of NGC 2403 the abundances of the five H II regions from Garnett et al. 
( 1997 ) at galactocentric distances larger than 2.7 kpc (0.25 R 2 s) were averaged. 



Abundances in NGC 4258 
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A.3. Oxygen abundances for Cepheid host galaxies using the transformation to the T e scale 

For galaxies in which no (or only a few) [O III] A4363 line detections are available we used Eq.[3]to scale the oxygen abundances 
calculated in the ZKH scale to the T e scale. The results are presented in Table [9] The galaxy list was drawn from Sak ai et aL] 



(20041, Saha et al. ( 2006[ > andjBono et al. (2010 1, with the addition of a few galaxies. Unless otherwise noted (see references at 
the end of the table), the abundances in the ZKH scale were taken from Ferrares e et al.| ( [2000l ). 



We comment below on a few galaxies from Table|9]in which auroral lines have been detected: 



NGC 224 (M31): A value of 12 + log(Q/H) = 8.41 ± 0.02 has been reported by Esteban et al. (2009 1 from a detection of 



[O III] A4363 in the outlying H II region K932 (R ~ 16 kpc from the galax y center). Given the rather shallow abundance gradient 
in this galaxy (still uncertain, but between -0.006 and -0.027 dexkpc~ 1 , |Trundle et al.|2002 i, this direct r c -based abundance is 
probably a good measurement of the overall metallicity. 

NGC 1365: |Bresolin et IE] ( |2003] > measured chemical abundances from the [N II] A5755, [S III] A6312 and [O II] A7325 auroral 
lines for three H II regions in the vicinities of the Cepheids studied by the HST Key Project. The weighted mean oxygen abundance 
is 12 + log(0/H) = 8.48 ±0.10, in good agreement with the value of 8.48 given in Table[9]for the T e scale. 

NGC 303 1 (M81): The detection of [Olll] A4363 in Hll regions in the galactocentric radial range 5.8-9.8 kpc has been re- 
ported bylStanghellini et al. (20101. The average abundance they derived is 12 + log(0/H) = 8.54^ 3|. This metallicity should be 
representative of the inner held in Table [9] which gives a lower value of 8.33. 

NGC 5236 (M83): |Bresolin et aT] pool] ) measured chemical abundances fro m the [N II] A5755 [S III] A6312 and [Oil] A7325 
auroral lines for two Hll regions in the vicinities of the Cepheids studied by |Thim et al.| p003). The weighted mean oxygen 
abundance is 12 + log(OZH) = 8.55 ±0.10, in good agreement with the value of 8.64 given in Table [9] 

NGC 7793: From the emission line data published by Edmunds & Pagel (1984) for two Hll regions with good [Olll] A4363 



detection we derive 12 + log(0/H) = 8.42 ± 0.09, in agreement with the 8.45 value in Table[9] 
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Table 9 

Metallicities of other Cepheid host galaxies in the T e and ZKH scales 
Galaxy 12 + log(Q/H) 





T e scale 


ZKH scale 




NGC55 


8.06 


8.35 a 


0.15 


NGC224 


8.49 


8.98 


0.15 


NGC925 


8.19 


8.55 


0.15 


NGC 1309 


8.44 


8.90 b 


0.19 


NGC 1326A 


8.16 


8.50 


0.15 


NGC 1365 


8.48 


8.96 


0.20 


NGC 1425 


8.50 


9.00 


0.15 


NGC 1637 


8.56 


9.08 c 


0.15 


NGC 2090 


8.37 


8.80 


0.15 


NGC 2541 


8.16 


8.50 


0.15 


NGC 2841 


8.71 


9.30 d 


0.15 


NGC 3021 


8.46 


8.94 b 


0.25 


NGC 3031 (inner) 1 


8.33 


8.75 


0.15 


NGC 3031 (outer) 2 


8.23 


8.60 e 


0.15 


NGC 3198 


8.23 


8.60 


0.15 


NGC 3319 


8.08 


8.38 


0.15 


NGC 3351 


8.67 


9.24 


0.20 


NGC 3368 


8.64 


9.20 


0.20 


NGC 3370 


8.37 


8.80 f 


0.19 


NGC 3621 


8.33 


8.75 


0.15 


NGC 3627 


8.68 


9.25 


0.15 


NGC 3982 


8.33 


8.75 f 


0.25 


NGC 4321 


8.59 


9.13 


0.20 


NGC 4414 


8.64 


9.20 


0.15 


NGC4496A 


8.35 


8.77 


0.15 


NGC 4527 


8.33 


8.758 


0.15 


NGC 4535 


8.64 


9.20 


0.15 


NGC 4536 


8.40 


8.85 


0.15 


NGC 4548 


8.74 


9.34 


0.15 


NGC 4639 


8.50 


9.00 


0.15 


NGC 4725 


8.45 


8.92 


0.15 


NGC 5236 


8.64 


9.2& 


0.15 


NGC 7331 


8.28 


8.67 


0.15 


NGC 7793 


8.45 


8.92 1 


0.12 


IC4182 


8.09 


8.40 


0.20 


References. — (a)|Zaritsky et al.|(1994|; (b)IRiess 


et al. 1 2009a i; (c) Leonard et al. 1 2003 i; (d) 



Saha et al. i 2006 1; (h) using line fluxes from Bresolin et al. i 2005 i; (i) using line fluxes from 




